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We present a static, model-independent, experimental determination of the spin-spin correlation 
length 5 in a quantum spin system. This is achieved in doped Haldane (i.e., 5=1 Heisenberg 
antiferromagnetic) chain Y2BaNii_i:Mgi:05 by *^Y NMR imaging of the staggered magnetization 
induced around the Mg-impurities (i.e., chain boundaries) by magnetic field. The magnitude of 
this magnetization is found to decay exponentially with ^ equal to the theoretical prediction for 
an infinite 5* = 1 chain, and the staggered magnetic moment at the edge site showing the Curie 
behavior of an effective S — 1/2 spin. 



PACS numbers; 75.10.Jm, 76.60.-k 

One dimensional (ID) quantum antiferromagnetic 
(AF) spin systems such as spin-ladders, spin-Pcicrls sys- 
tems and spin chains have recently attracted an increas- 
ing interest, partly as a possible way to approach from a 
lower dimension the physics of the copper-oxide high Tc 
superconductors. In all these systems the spin-spin corre- 
lation function is a quantity of fundamental importance, 
which can be experimentally probed as a function of wave 
vector and energy by inelastic neutron scattering, or by 
studying the static response of the system to the pertur- 
bation created by impurities, through NMR imaging of 
local spin polarization. In recent examples of this latter 
technique applied to high Tc superconductors ||l|, spin- 
ladder systems g] or 5 = 1/2 chains Q, the broadening 
of the NMR lines (or the modification of the lineshape) 
due to impurities has always been interpreted within a 
specific model in order to obtain information on the cor- 
relation length ^, making conclusions manifestly model 
dependent. In this letter we present an NMR study of a 
doped Haldane (i.e., 5 = 1 Heisenberg AF spin) chain, 
in which individual spin sites near impurities are resolved 
as satellite peaks in the NMR spectra. In this way the 
real-space dependence of the spin polarization is followed 
site by site, providing a complete and model independent 
picture of the impurity screening, and giving thus a di- 
rect access to ^. Such a study should be compared to the 
observation of satellite lines revealing screening of normal 
M or Kondo m impurities in a metal, but it applies here 
to a low-dimensional quantum antiferromagnet, i.e., to a 
class of systems in which screening of non-magnetic im- 
purities has generated recently a large body of theoretical 
works Q. 

The main characteristic of integer spin chains is the 
existence of a gap between the lowest energy levels, pre- 
dicted by Haldane |7|| and experimentally observed in sev- 
eral compounds ||q-]l^ . If such a system is doped by non- 



magnetic impurities, an ideal infinite chain is "cut" to an 
ensemble of open segments of length L, in which the low- 
est lying eigenstates are a singlet |0) and a triplet |1, S^), 
with an energy separation exponentially decreasing with 
L Q. The Haldane gap, Q defined for L -^ oo, is the 
one above the resulting fourfold quasi-degenerate ground 
state. In the state |1, 1) it has been shown by numerical 
calculations [p5l-p8l , that the expectation value of the lo- 



cal spin operator (S^) decreases exponentially 19 1 from 
the boundary with staggered phase. In particular, for a 
pure Heisenberg model, one has an effective spin-1/2 at 
the edge of the chain_((S'f ) ~ 1/2) and a characteristic 
decay length ^ ~ 6 |1^, a value very close to the zero- 
temperature correlation length (^oo) of an infinite chain. 
Thus, it has been argued that the zero-temperature local 
magnetization (corresponding for long chains to (S^) in 
the state |1, 1)) reflects indeed the decay of the correla- 
tion function in the Haldane phase. 

S* = 1/2 boundary spins have been experimentally 
observed in low-temperature ESR spectra of pure and 
doped Ni(C2H8N2)2N02(C104) |^,|l| and YaBaNiOs 
[ p2[ . Also specific heat measurements psf c an be ex- 
plained in terms of 5 = 1/2 excitations ]24[|. On the 
other hand, no direct experimental determination of the 
site dependence of the local magnetization have been re- 
ported until now. Moreover, the temperature dependence 
of the expectation value of edge spins, and the relation- 
ship between the decay of the staggered magnetization 
and the infinite chain correlation function have not been 
established yet at temperatures where the excited states 
are populated. 

The compound Y2BaNi05, due to the low in- 
ter/intrachain coupling ratio {\J'/J\ < 5 x 10""*, with 
J = 285 K) and the small single- ion anisotropy [pdj , 
is a prototype one-dimensional 5 = 1 Heisenberg AF. 
Moreover, a controlled substitution of Ni^+(5 = 1) by 



Mg^+ {S = 0) offers the possibility to obtain open chains 
of the average length (L) = 1/x. In Y2BaNii_a:Mga;05, 
Ni^+ local magnetization generates a site dependent hy- 
perfine field at *^Y nuclei inducing a multi-peak struc- 
ture in *^Y NMR spectra. These spectra are sort of a 
map of (Sf), inaccessible by any other technique because 
a non-homogeneous, variable- amplitude response to the 
external field is involved. 

*^Y NMR measurements have been carried out on a; = 
0.05 and x = 0.10 powder samples of Y2BaNii_j:Mga;05, 
at the NMR facility of the Grenoble High Magnetic Field 
Laboratory. The ^^Y NMR spin-echo spectra have been 
obtained at fixed frequency by sweeping the magnetic 
field. Fig. 1 shows representative spectra for the x — 0.05 
sample. Several peaks can be observed, their shift from 
the central line increasing with decreasing temperature. 
In the following we show that this temperature depen- 
dence yields directly the temperature evolution of the lo- 
cal magnetization. The shift A(p/) = \Ho{pi) — Ho{pc}\ 
of the edge-spin peak pi=i with respect to the central 
peak Pc is reported in Fig. 2 as a function of T. A 
comparison of spectra for two samples with different av- 
erage chain lengths shows that for T ~ J the peak po- 
sitions are independent on (L) (see the inset to Fig. 2). 
Fig. 3a shows the position dependence A(p/) at differ- 
ent temperatures. It is found that the shift decreases 
exponentially with / at all temperatures above ~ 100 K. 
Below 100 K the peaks are smeared into a single broad 
line (inset), preventing us to extend the analysis to lower 
temperature. 

The intensity of the NMR spectrum is proportional 
to the number of ^^Y nuclei which, for the irradiation 
frequency vrf, obey the resonance condition 
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where (7/277) = 2.0859M Hz/Tesla is the ^^Y gyromag- 
netic ratio, M the z-component of the hyperfine field at 

the j-th nucleus and Aikj the Ni^+-®^Y hyperfine cou- 
pling tensor. In Eq. (1) fc is a chain index and / is 
an in-chain site index. A given ^^Y has two Ni^+ near- 
est neighbors (nn) belonging to the same chain and two 
next nearest neighbors (nnn) on adjacent chains ||2^. On 
the basis of the NMR results in pure Y2BaNi05 one can 
assume that only nnn sites contribute to h^ : 

A"/^ ^ A = 13 kG ii a — P and (j) - (/, fc) are nnn 



^ikj 
At = 



if a 7^ /3 or (j) - {I, k) are not nnn 

(2) 



tra jl^ having a shift much larger than the line broad- 
ening indicate that the hyperfine coupling is essentially 
scalar and hence due to a transferred hyperfine interac- 
tion. Since the mechanism of transferred hyperfine inter- 
action requires extension of the Ni^+ — 3d orbitals at the 
Y site, only nn and nnn Ni^+can couple to ^^Y. On the 
other hand, as already noticed |12], the lattice position of 
the *^Y is strongly unfavorable for transferred interaction 
from nn Ni^+, and the coupling can be assumed essen- 
tially limited to two equivalent contributions from nnn 
ions only. Since for an infinite chain in uniform magnetic 
field the magnetization is uniform, {Sik) ~ {S), from the 
standard plot of the shift of the NMR line as a function 
of the spin susceptibility one derives J^i k A-ikj = 2A = 
26 kG. Furthermore, small interchain coupling J' allows 
us to consider the nnn ions (which belong to different 
chains) contributing to M as completely uncorrelated. 
Thus, Eqs. (1) and (2) link the position of the max- 
ima in the NMR spectra to the different values of the 
local magnetization along one chain. In particular, the 
high temperature spectra shown in Fig. 1 correspond to 
a large number of Ni^"*" spins having the same expecta- 
tion value (Sf) = {S^) (central peakpc), and to a certain 
number of sites for which (S^) ^ {S^), yielding the peaks 
Pi . According to T = numerical results |15|-[l^] , one ex- 
pects \{Si) - {S^)\ -> for i ^ L/2. Thus we map the 
peaks Pi with the positions / and (L — / + 1) along the 
chain, e.g., positions 1 and L correspond to pi. 

As shown in Fig. 2 (solid line) the temperature evolu- 
tion of A(pi) is very well accounted for by a hyperfine 
field of the form 



[A{p,)]iT)^A 
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In the following we summarize the main reasons sup- 
porting this assumption [p6|. The NMR powder spec- 



with S = 1/2, g = 2 and Hq = 14 Tesla. This means 
that, for A/f < T < J (where Ah « 100 K is the Haldane 
gap for pure Y2BaNi05), the non-uniform magnetization 
associated to the extreme sites of the chain is close to 
the one corresponding to "free" S = 1/2 spins. Since no 
doping (i.e., (L)-) dependence of A(pi) have been found 
at room temperature (inset to Fig. 2), this conclusion 
is independent on the chain-length at least for (L) > 10 
and T ~ J. One could observe that a, S = 1/2 Curie-like 
behavior can be obtained as thermal expectation value 
of the edge spins on the lowest energy states (singlet and 
triplet) of the finite chain, provided they are separated 
by an energy gap AE << gfisHo. According to Monte 
Carlo calculations |1^ , AE is indeed considerably smaller 
than the Zeeman energy in a field of 14 T (~ 0.08./) 
for L ~ 20 {i.e., x = 0.05), while AE « .g/isiJo for 
L ~ 10 {i.e., X = 0.10). However, one has to notice that 
for T > Ah a precise calculation of the thermodynamic 
properties of the system has to take into consideration 
also the excited states, as is evident from the analysis of 
the spin-spin correlation reported below. 



The dependence A(p/) on / can be fitted to an expo- 
nential decay 



[A{pj)]{T) = [A{p,)]{T)< 



-(7-l)/4(T) 



(4) 



(solid lines in Fig. 3a), and the characteristic decay 
length £,{T) of the magnetization can been extracted. As 
shown in Fig. 3b, ^(T) obtained in this way is very close 
to the values of the infinite chain correlation length ^oo 
calculated by Kim et al. p^. This remarkable conclusion 
indicates that the picture of boundary defects, which re- 
flects the bulk correlation of an infinite chain, holds not 
only in the Haldane phase [T « Ah) fl^], but in gen- 
eral for T < J. A possible way to interpret this result 
can be based on a heuristic extension from zero to finite 
temperature of the model proposed in Ref. [|7|- There 
the expectation values of Sf in the state |1, 1) are shown 
to decay exponentially over a length ^ = S^^(T — 0), due 
to the interaction of the fields describing the translational 
invariant chain in the free boson model, with the effective 
fields originating from two end of chain S — 1/2 spins. 
The experimental findings reported in Fig. 3 are then im- 
mediately recovered if one assumes that these edge " ex- 
tra" spins are practically Curie-like (as indicated by the 
behavior of pi in Fig. 2), and substitutes £,oo{T — 0) by 
£,oo(T) and the expectation values of 5*^ in |1, 1) by the 
thermal average. 

At low temperature the progressive increase of the 
width d{pi) of the NMR peaks (see Fig. 1) prevents a 
meaningful analysis of the boundary magnetization along 
the lines used for T > 100 K. However, the characteristic 
shape of the NMR spectra pq| clearly indicates the per- 
sistence of field-induced staggered magnetization. The 
origin of the broadening is likely related to the anisotropy 
terms in the magnetic Hamiltonian nW , which in powder 
samples causes a distribution of (Sf). This observation 
is corroborated by the fact that the ratio between the 
easy-axis anisotropy D = 8 K and the Zeeman energy is 
of the same order of magnitude as d{pi)/A{pi). In order 
to extend the data to lower temperature it is thus neces- 
sary to work on single crystals. However, an additional 
source of broadening is a distribution of chain lengths 
around the average value (L), which is expected to be 
important at T < 100 K. Indeed, at low temperature ^ 
is increasing, and for short chains with L/2 < ^ mag- 
netization is becoming L-dependent. In principle these 
effects can be resolved by decreasing the dopant concen- 
tration X, however, this will proportionally decrease the 
size of our satellite peaks, i.e., the "resolution" of NMR 
spectrum. 

In summary, through *^Y NMR in Mg doped 
Y2BaNi05, we have been able to image the local mag- 
netization induced by "free" boundaries in a S* = 1 
Heisenberg AF spin chain, in the temperature range 
Ah ^ T '^ J. By mapping the lattice positions onto 
the different peaks in the NMR spectra, we have shown 



that the characteristic decay length ^(T) of the bound- 
ary magnetization is equal to the theoretical prediction 
for the correlation length of the Haldane chain. The non- 
uniform magnetization of the edge spin in a finite chain 
follows a 5* = 1/2 Curie-like behavior. These results pro- 
vide direct experimental confirmation that the NMR of 
doped quantum AF spin systems is the right technique to 
access the spin-spin correlation function. Improving the 
choice of experimental system {e.g., by using single crys- 
tals) and using an appropriate nucleus for NMR imag- 
ing, this technique should provide more reference data 
for comparison to the theory. 
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FIG. 1. *^YNMR spectra in Y2BaNio. 95 Mgo.osOs recorded 
at fixed frequency urf = 29.4 MHz by sweeping magnetic 
field. Resolved satellite peaks are labelled with the index J, 
following the decreasing magnitude of their shift (measured 
from the central line). In the inset, all the peaks are shown 
to be smeared in a single wide line when the temperature is 
lowered. 



FIG. 2. Temperature dependence of the shift of the pi peak 
in *^Y NMR spectra (measured from the central line, see Fig. 
1). The solid line is a Curie function for S = 1/2 spins (Eq. 
4 in the text). Inset shows the doping dependence of the 
room temperature spectra of Y2BaNii_2,Mg2,05 demonstrat- 
ing that peak positions do not depend on the chain length. 



FIG. 3. a) Site-dependence of the magnitude of shift of 
Pi peaks in Y2BaNio.95Mgo.05O5 taken at prf = 29.4 MHz. 
Solid lines are fits (Eq. 5 in the text), defining a decay length 
5(r) of the boundary magentization. b) Comparison of an 
experimental ^(r) value (full squares) and the infinite chain 
correlation length ^oo{T) (up triangles) as calculated in Ref. 
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